Interlayer couplings and the coexistence of antiferromagnetic and d-wave pairing 

order in multilayer cuprates 
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A more extended low density region of coexisting uniform antiferromagnetism and d-wave super- 
conductivity has been reported in multilayer cuprates, when compared to single or bilayer cuprates. 
This coexistence could be due to the enhanced screening of random potential modulations in in- 
ner layers or to the interlayer Heisenberg and Josephson couplings. A theoretical analysis using a 
renormalized mean field theory, favors the former explanation. The potential for an improved de- 
termination of the antiferromagnetic and superconducting order parameters in an ideal single layer 
from zero field NMR and infrared Josephson plasma resonances in multilayer cuprates is discussed. 



I. INTRODUCTION 



Although the cuprate superconductors have been stud- 
ied for a quarter century, the exact form of the interplay 
of antiferromagnetic (AF) and d-wave pairing(dSC) order 
at low hole doping is not known reliably. Experiments on 
very low doped cuprates show insulating behavior and a 
critical hole concentration, Xc, for the onset of supercon- 
ductivity. The presence of a strongly varying potential 
due to the random distributions of the acceptors is be- 
lieved to be important at very low hole densities. There 
are reports of a spin glass phase separating the AF and 
dSC regions of the phase diagram^. A series of neutron 
scattering experiments on La2-a;Srj;Cu04 found stripe 
order but with a change of orientation from parallel to the 
Cu-O-Cu bonds at a; > Xs to a 45° direction at a; < x^. 
Neutron scattering experiments on YBa2Cu306-i-j/ found 
spin density wave (SDW) order in the hole concentration 
range 0.05< x < 0.08'^. On the other hand coexistence of 
uniform AF order with dSC has been deduced from NMR 
studies of multilayer Hg- and Ba2Ca3Cu408(FyOi_y)2- 
cuprates in a substantially larger hole density ranged . In 
this paper we analyse the role of AF and dSC interlayer 
coupling in multilayer cuprates and their effect on the 
coexistence on AF and dSC order. 

The observation by Kitaoka, Mukuda and 
collaborators^, of substantial zero magnetic field shifts 
on a unique Cu-site in the innermost layers led them 
to claim uniform AF order in underdoped layers with 
X ^ 0.1. The multilayer cuprates are doped by acceptors 
in the insulating blocks between the multilayer blocks, 
but their random potential will be weakened in the 
inner layers due to screening by the more highly doped 
metallic outer layers. Thus these multilayer cuprates 
offer the best possibility to reliably determine AF and 
dSC ordering at low doping in a clean single layer. 
In this paper we illustrate the potential of multilayer 
systems by considering the case of 4-layer cuprates 



for which experimental results exist in a wide density 
range. Variational Monte Carlo calculations for an ideal 
single layer strong coupling t-J model find a substantial 
coexistence region of spatially uniform AF and dSC 
order up to a critical hole density x ^ 0.1^-. Theoretical 
support for dSC order at low doping in a single layer 
also comes from exact diagonalization studies on clusters 
of the t-J model containing up to 32-sites doped with 
2 holes^. When extrapolated to an infinite layer, these 
point to a robust d-wave cooperon ( bound hole pair 
) resonance at low doping . The virtual exchange of 
cooperons will act as a pairing mechanism for doped 
holes in near nodal states at low hole densities, in the 
presence of strong AF local (and probably long range) 
order^^. 

The two order parameters to be determined are the 
sublattice magnetization of AF order and the pairing am- 
plitude of dSC ordering. The former is directly measured 
by the zero field shift on Cu-sites in NMR experiments 
In multilayer samples the Josephson couplings between 
the dSC order in neighboring layers lead to the Josephson 
plasma resonances which appear in infrared spectra with 
the electric field oriented along the c-axis'^. The energies 
of these resonances are proportional to the product of 
the dSC order parameters. The combination of these two 
measurements on the same samples can, in principle, be 
used to determine the values and the dependence of both 
AF and dSC order parameters on the hole density. To 
obtain information on the interplay of AF and dSC order 
in a single layer, we need to include the effects of inter- 
layer couplings between the respective order parameters 
in the analysis. Note, the interlayer spin coupling dis- 
favors SDW order in adjacent planes with different hole 
densities and therefore different periodicities. As a result 
we restrict our attention to the case of commensurate AF 
and dSC order. 

In this paper we report on calculations using the renor- 
malized mean field theory (RMFT), a method introduced 
in the early days by Zhang and collaborators^, to treat 
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the strong coupling t-J model. We extend the method 
to multilayer t-J models. RMFT calculationsi2 generally 
agree well with more accurate variational Monte- Carlo 
calculation»ii in the case of a single layer. 

To treat multilayer materials we require the values 
of the interlayer couplings. The Heisenberg spin-spin 
coupling has been directly measured with a strength 
typically of order a tenth of the intralayer coupling. 
The interlayer Josephson coupling constants are less well 
known, as will be discussed further below. 

A second key input is the value of the hole densities 
in inequivalent layers. Kitaoka and collaborators have 
measured the Cu-Knight shift in each layer and used the 
room temperature value as input in an empirical formula 
to estimate the hole density in each layer. We shall return 
to these inputs below. 
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with uJc a characteristic energy. This excited state 
Cooper pair can subsequently relax to the Fermi sur- 
face of the second layer through interactions such as the 
J-term in the single plane Hamiltonian. This virtual tun- 
neling process through pair states off the Fermi energy 
will be the leading contribution of the interlayer Joseph- 
son coupling, arising from Cooper pair hopping between 
inequivalent layers leading to an interlayer Josephson 
coupling 



II. t-J MODEL FOR MULTILAYER CUPRATES 



We start by defining the t-J model for a single plane. 

Tt-J 



Hi for the I layer is given by 
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where cua- is an annihilation operator of a spin a elec- 
tron at site i in layer I, Pq is the Gutzwiller projection 
operator which enforces the no-double occupancy condi- 
tion in cuprates, tij is the hopping integral between site 
i and site j with nearest neighbor hopping tjvA' — t, next 
nearest neighbor hopping t at at at = t\ and tij = for 
the remaining site pairs. J is the superexchange coupling 
between NN sites. 

The interlayer coupling, iJim, between two neighbor- 
ing layers I and m includes a superexchange spin-spin 
coupling and an interlayer hopping term. The interlayer 
superexchange coupling is given by 

Hiin ^ J^'Y • Sm,i, (2) 

i 

and the interlayer hopping term reads 

HL =PgY. i±(k)(c[k.c™,ka + h.c.)PG, (3) 



cos ky. 



wherein — t±(j)1^/2, with (/)k = cosfc^: 

In a system with more than two layers, the hole concen- 
trations on two adjacent inequivalent layers are different 
and the Fermi surfaces of the two layers are mismatched. 
Because planar momentum is conserved in interlayer hop- 
ping, the mismatch of the Fermi surfaces means that di- 
rect interlayer single particle hopping and interlayer pair- 
ing can be neglected. However a Cooper pair can tunnel 
from the Fermi surface of one layer into a pair state off 
the Fermi surface on a neighboring layer, conserving the 
planar momentum of the individual electrons. This pro- 
cess can be written down in standard perturbation theory 
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X Pg (^c[-ki4kt'^™,k'tCm -k'i + h.c.j Pg, (4) 

where J(q) = J{cosqx + cosqy) is the pair hopping am- 
plitudes, uj'^ is a characteristic energy related to the in- 
tralayer relaxation process, and the prime on the summa- 
tion denote s that the summation are only in the vicinity 
of Fermi surfaces with cut-off ujd ■ 
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FIG. 1: Schematic illustration for the hole carrier charge dis- 
tribution X, electric field E (direction indicated by arrows) 
in 4-layer Cu-oxides. Horizontal lines represent layers. Sub- 
indices i and o represent inner and outer planes, respectively. 
The top and bottom dashed lines represent chemically doped 
charge. 

In multilayer materials the distribution of hole den- 
sities on the individual layers plays an important role. 
Theoretical estimates of the hole densities can be made 
using a Hartree approximation scheme by including the 
electrostatic energy in the Hamiltoniani^. For simplicity, 
we consider each layer as an infinite plane with homoge- 
neous charge. The electric field generated by such a plane 
75— ^^T, where is the relevant dielectric 



is £; = ^ = 



constant, x is the hole concentration of the layer, and a is 
the lattice constant. For a 4-layer system, we denote Xo 
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and Xi as the hole concentrations of the outer layers and 
the inner layers, respectively. The average hole concen- 
tration is a; = {xo + Xi)/2, while the electron concentra- 
tion of the acceptor layer above and below the 4-layers is 
2x. The spatial distribution of electric field is depicted 
in fig. [TJ We are only interested in the terms relevant 
to the charge inhomogeneity, so the electrostatic energy 
between acceptor layer and outer layer can be ignored. 
Then the electrostatic energy per plaquette in the 4-layer 

case is -Ep, = where d is the distance between two 

adjacent Cu02 layers. The Hamiltonian for multilayer 
system then reads 



I {Irn) 



Him) + 



(5) 



The dielectric screening of the interlayer electric field 
arises from displacements of the ionic layer between the 
neighboring layers and is not identical to the static c- 
axis dielectric screening due to optical phonons etc. We 
adjust the value of to fit the estimated values for the 
layer densities (see below). 



We now introduce the following mean fields 
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where ()g indicates average with the unprojected state. 
Since we are only interested in the ground state, 
we can assume that the mean fields Aq,/ are all 
real. In the presence of AF order, the magnetiza- 
tion on the sublattices A and B are opposite, and 
the mean fields XQ.i,ij,a from NNN sites depend on 
the sublattice and spin. We introduce the following 
mean fields, xo,i.AA.a = jNT,teAT,j=NNN{z)Xo,i.ij.a 
and Xo,i,BB,a = -mT.ieBJ2j=NNN(i)Xo..i,ij,a- Since 
Xo,LAA,a = Xo,i,BB~a, we Can further simplify the no- 
tation as Xo,i,AA = Xo,i,AA,t and Xo,i,bb = Xo,LBB,t- On 
the other hand, the mean fields XQ,i,ij.(7 for NN sites are 
independent of the sublattice ij and spin cr, so we can 
set the mean fields Xo,i,ab = m X^^a T.j=NN{i) Xo,i,ij,a- 
The hole concentration on each layer satisfies the condi- 
tion. 
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III. RENORMALIZED MEAN FIELD THEORY 
FOR MULTILAYER SYSTEM 



In this section, we will use RMFT to analyse the multi- 
layer Hamiltonian. The main element of the RMFT the- 
ory is the so-called Gutzwiller approximation, which re- 
places the Gutzwiller projection in the Hamiltonian by 

simple numerical factors. Let and be the av- 

erage values of operator O in the Gutzwiller projected 
state and in the unprojected state, respectively, then 

~ 5o 1 '^ith go a numerical factor, depend- 
ing on the hole density and the process associated with 
the operator O. The coexistence of AF and dSC order in 
the t-t'-J model for a single layer material has been in- 
vestigated within RMFT recently by K.-Y. Yang et alJ^. 
Here we generalize their calculations to a multilayer sys- 
tem. Note that in the RMFT the superconducting order 
parameter is non-zero at wave-vector off the Fermi sur- 
face. The virtual tunneling process through pair states 
off the Fermi energy when the two layers have different 
hole concentrations we discussed in Section II may also be 
taken into account by including the pair scattering on the 
same layer along with the second order interlayer hopping 
process conserving momentum of single electrons. Below 
we shall use Eq. (4') for the interlayer tunneling directly 
to study the multi-layer system described in Eq. (5). 



Then we consider the Gutzwiller factors. In in- 
tralayer terms, we adopt the same form of the factors 
9f,AB^9*i,AA^9f,BB:9r:9r^ and gf as used previously for 
the single layer t-J model^". For the interlayer hop- 
ping and superexchange coupling terms, we make a sim- 
ple assumption that the Gutzwiller factors are g*,^ = 



91ab9Lab and gf. 



\j9T9nn respectively. The to- 



tal energy reads 



E — E + E ^ 

I ilra) 



(8) 



where £"/ is given by the RMFT for single layer t-t'-J 
modeli^, the electrostatic energy Ef,g is specified in the 
previous section, and the interlayer energy Ei^a takes the 
form 

Elm = ^"^'^ + = _4(.9L)'T^E4(Ao.,kAo,™,k 



Ao,;,k+QAo,m,k+Q) + J±9lm^0,l'mn 



(9) 



with Ao,;,k = y-kfC-i^ij ■ Note that the magnetic mo- 
ment mi and superconducting order parameter A/ are 
related to the mean fields by the Gutzwiller renormaliza- 
tion factors^. 
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^i^91ab^o,i- (10) 
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These approximations lead to the mean field Hamilto- 
nian 



- X! (^,kC/,-k4,cz,kt + h.c.) , (11) 



No- 
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and 7k = cos + cos ky , and 0k — cos fc^; cos fcy . 



IV. RESULTS AND DISCUSSIONS 
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FIG. 2: The order parameters and hole concentrations calcu- 
lated in 4-layer t-t'-J model Eq. (5). Model parameters are 
J = 0.3t, t' — ~0.3t, interlayer pair hopping t^/cjc = O.lt, 
interlayer spin coupling J± = O.IJ. Upper panel: the hole 
concentrations for each layer xi vs average hole concentra- 
tions X. Thick blue curves are for inner layer (xi), thin red 
for outer layer (xo)- Lower panel: magnetic moment m and 
superconducting order parameter A as function of hole doping 
for each layers xi. Thick blue curves are for inner layer, thin 
red for outer layer. Left panel: dielectric constant e = 200, 
right panel: e — 50. Dashed curves in (b) and (d) are for 
single layer results, plotted for comparison. 



The mean field energy defined by Eqn (8) has been 
minimized numerically for the case of 4-layers. The re- 
sults depend on the hole densities which in turn depend 
on the value of the dielectric constant, e^. In fig. [H the 
two upper figures show the density in each layer as func- 
tion of the total doping for two choices — 200 and 
er = 50. The lower figures show the resulting values of 



the order parameters in each layer. The single layer result 
is also shown in the figure as a dashed line for compari- 
son. The SC order parameters of the 4-layer system and 
single layer system are almost equal and independent of 
er ■ This is a consequence of the weak Josephson coupling 
between adjacent layers in the multilayer system due to 
the mismatch of the Fermi surfaces and the suppression 
of the interlayer hopping by the strong onsite Coulomb 
interaction. 



TABLE I: Comparison between NMR experimental results on 
Ba2Ca3Cu408(Fj;Oi_y)2^ and the renormalized mean field 
theory for 4-layer t — t' 
in the compound, x 
on the outer and inner planes, m 
magnetic moment on the outer and inner planes, respectively. 
e is the fitted value of the dielectric constant from the hole 
distribution, lju and ujio are the Josephson plasma frequencies 
between two inner layers and between inner layer and outer 
layer, respectively. 



— J model, y is the content of F 
op and Xip are the hole concentrations 
op and rriip are the AF 
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The interlayer AF coupling slightly enhances the AFM 
moment and critical doping Xc of the inner layers and 
leads to a finite AFM moment on the outer layers when 
their hole concentration is much larger than the single 
layer critical doping Xc- Such a long tail arises because 
of the lower hole concentration on the inner layers. A 
smaller leads to larger charge imbalance between the 
outer and inner layers and so to larger Xc of outer layers, 
as shown in fig. [5] If one neglects the long tail of the 
AFM moment of the outer layers, the phase diagram of 
multi-layer system is very similar with that of a single 
layer system, independent of the value of the dielectric 
constant e^- This indicates that experiments on multi- 
layer system can capture the essential physics of single 
layer t-J model 

We turn now to the relation of the model with ex- 
periment. The two order parameters we investigated 
above are the AFM magnetization and the pairing am- 
plitude of dSC ordering. The former has been directly 
measured by the zero field shift in NMR. In the NMR 
experiments, Shimizu, Kitaoka and collaboratorsi^ have 
also used an empirical scaling of the Cu Knight shift, 
measured at room temperature, to deduce hole densities 
of individual layers. Note this scaling form has evolved 
with time. In tab. HI we show the value of the rela- 
tive dielectric needed to fit the charge distribution re- 
cently determined from NMR experiments on the 4-layer 
Ba2Ca3Cu408(Fj^Oi_y)2 cuprate. Note, in these oxy- 
fluoride cuprates the hole density differs substantially 
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from the nominal F- concentration, y. Though the dielec- 
tric constant required to fit the hole density changes from 
87 to 247 as the total hole density decreases from 0.46 to 
0.264, the physics in each layer should not change very 
much as analyzed above. In Table 1, wc also show the 
AFM moments calculated by the RMFT and measured 
by NMR experiments. The former values are consider- 
ably larger than the latter, indicating an over estimation 
of the AFM order in the RMFT theory. The comparison 
between theory and experiment is sensitive to the values 
of the hole density. For example, if we scale the value 
of the estimated hole density so that the optimal density 
with the highest Tc for a single layer is {xc = 0.2) rather 
{xc ~ 0.16) quoted in Mukuda et al^^, this leads to an 
increase in the input values for the hole density in the 
RMFT calculations by a factor of 1.25. Then the cal- 
culated values of the AF moment shown in Tab 1, are 
smaller, although still too big, and in better agreement 
with experiment. 

Turning to the dSC ordering, as mentioned earlier the 
pairing amplitudes enter into the Josephson plasma fre- 
quencies observed in the infrared spectra with the electric 
field oriented along the c-axis. In multilayer cuprates, 
there are two kinds of Josephson plasma frequencies 
which correspond to inter-multilayer and intra-multilayer 
Josephson couplings respectively. In this paper, we are 
only interested the latter, i.e. the Josephson plasma fre- 
quencies related to the Josephson coupling between two 
adjacent Cu02 layers in same unit cell. The Joseph- 
son energy between those two layers can be written as 
Ej = — Jim cos{(f>i — 4>m), where 4>i and 4>m are the phase 
of the superconducting order parameter of the two lay- 
ers respectively. For the two layers with the same hole 
concentration, J;,„ is given in our RMFT by. 



(13) 



where J-i is the anomalous Green's function of layer /, 
iijJn is the Matsubara frequency, Q = (tt, tt), and indi- 
cates the summation is only over the reduced zone. The 
Josephson plasma frequency which corresponds to the 
Josephson coupling between those two layers is propor- 
tional to the Josephson coupling constant^ii^, i.e. 



So we have 



(14) 



where C is a constant related to t_L and finite frequency 
dielectric constant e(wj). For simplicity, we can assume 
that C is weakly dependent on the hole concentrations 
and materials. 

For a 4-layer material, there are two intra-unit cell 
Josephson plasmon frequencies ojio and uja which cor- 



respond to the Josephson coupling between two inner 
layers and between the inner and outer layers respec- 
tively. A recent optical experiment on a 4-layer Hg-based 
compound has shown that the Josephson plasmon fre- 
quencies are 360 cm~^ for inner layers and 540 cm^^ for 
outer and inner layers respectively leading to a ratio of 
1.5^^. According to eqn. (IT4l) . ujiojion = oJio/oJu- Using 
the estimated hole concentrations of the Hg-compound, 
the ratio calculated with RMFT is around 1.25 which 
is fairly close to the experimental value. To make fur- 
ther comparison, we calculate the factor C with C — 
^JujioUJii/CjiobJii « 10, where the plasmon frequencies are 
the experimental values for the Hg-compound and the u 
is calculated with RMFT. Using this estimate for C al- 
lows us to obtain values for the plasma frequencies for 
the 4-layer Ba2Ca3Cu408(FyOi_y)2 material. The re- 
sults are quoted in Table HI 



V. SUMMARY 

In this paper we presented a series of calculations using 
the RMFT to treat the AF and dSC order parameters in 
4-layer cuprates. The RMFT calculations/citestripe-rvb 
and also the results of single VMC calculations/citeVMC, 
give a more extended region of coexistence of the two 
order parameters than that observed in single and bi- 
layer cuprates. The extended coexistence region in the 
4- layer cuprates agrees with the conclusion of the NMR 
experiments in the inner layers. In our analysis the ex- 
tended region comes not from strong intcrlayer magnetic 
and Josephson coupling, but from the extended coexis- 
tence region already present in single layers in the RMFT 
approximation. This analysis suggests that the differ- 
ence in the density range of coexisting order in multi- 
layer cuprates compared to experiments on single and 
bilayer cuprates, is due to the better screening of the 
external potential modulation in the inner layers of mul- 
tilayer cuprates, when compared to single and bilayer 
cuprates. The latter are adjacent to the random accep- 
tors, while the metallic outer layers screen the random 
potential at the inner layers. In principle the combina- 
tion of zero field NMR and infrared Josephson plasma 
frequencies allows one to determine both AF and dSC 
order parameters. However, in order to go further and 
directly compare measured and calculated order param- 
eters a better determination of the interlayer Josephson 
coupling constant and its variation with the hole densi- 
ties in the individual layers, is required. To this end a 
series of measurements of the Josephson plasma frequen- 
cies and the AF moments by zero field NMR at different 
hole densities on the same multilayer material would be 
helpful. 
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